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Abstract. In this paper we present a unified approach to the spectral analysis 
of an hypergeometric type operator whose eigenfunctions include the classical 
orthogonal polynomials. We write the eigenfunctions of this operator by means 
of a new Taylor formula for operators of Askey- Wilson type. This gives rise to 
some expressions for the eigenfunctions, which are unknown in such a general 
setting. Our methods also give a general Rodrigues formula from which several 
well known formulas of Rodrigues type can be obtained directly. Moreover, 
other new Rodrigues type formulas come out when seeking for regular solutions 
of the associated functional equations. The main difference here is that, in 
contrast with the formulas appearing in the literature, we get non-ramified 
solutions which are useful for applications in combinatorics. Another fact, 
that becomes clear in this paper, is the role played by the theory of elliptic 
functions in the connection between ramified and non-ramified solutions. 

Introduction 

A large number of discrete problems in graph theory, group representation theory, 
discrete harmonic analysis or finite symmetric spaces involve the spectral analysis of 
certain difference operators. Namely, given a positive integer d and denoting by 
the space of functions / : {0, 1, . . . , d} — ► C, we are interested in the eigenfunctions 
of linear operators A : V<j — > of the following type 

A/(f) = a+(t)(f(t + 1) - f(t)) - *_(*)(/(*) ~ /(* " 1)) 

where <j± £ V<j and cr_(0) = 0+(d) = 0. For instance, the work illustrates how 
this kind of eigenvalue problems become useful to solve some a priori unrelated 
combinatorial problems. Let us write ¥k[x] for the space of complex polynomials of 
degree < k in one variable. Explicit solutions of this spectral analysis are available 
when there exists an injective mapping x : t € {0, 1, . . . , d} t— * x% £ C such that the 
following conditions hold 

(a) We have Xt+i + %t-i = f(%t) for some / £ Pi [a:] and 1 < t < d — 1. 

(b) The subspaces Vd,k — {f °x : f € Pfc[s]} are A-invariant for < k < d— 1. 
By a well-known result of Leonard [5] -see also the reference [S]- we have (b a) 
whenever d > 7. When this is the case, the eigenfunctions of A can be obtained by 
applying the theory of basic hypergeometric series. Alternatively, Nikiforov, Suslov 
and Uvarov study in jTl] the spectral decomposition of A using properties (a) and 
(b) and some specific methods. This gives rise to a whole new approach to the 
theory of classical orthogonal polynomials. Following any of these methods, two 
very relevant expressions for the eigenfunctions of A can be obtained. First, as we 
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have recalled above, we can write the eigenf unctions in terms of g-hypergeometric 
polynomials. Second, the eigenfunctions can be written in general by means of a 
Rodrigues type formula. If we look at the approach which uses the theory of basic 
hypergeometric series, a long list of subcases must be checked. For instance, the 
reader is referred to the book of Bannai and Ito ^ . Besides, in that reference the 
authors do not provide any Rodrigues type formula for a problem which in essence 
is equivalent to the one considered here. On the other hand, in the approach of 
|14|. the derivation of Rodrigues type formulas is cumbersome. Moreover, at the 
very final stages, the authors need to invoke the Watson-Sears transformation and 
other formulas of the theory of g-series to get the hypergeometric expressions. 

In this paper we propose an intermediate approach which, from our point of 
view, is quite simpler and provides a somehow natural classification of the whole 
problem which seems to be unknown. Throughout the paper, we shall also obtain 
some results and remarks which, as far as we know, are new in the theory. Following 
Nikiforov, Suslov and Uvarov, we shall study the eigenfunctions of A with the aid 
of properties (a) and (b), so that the theory of basic hypergeometric series does 
not play a central role. However, in contrast with |14| . we reformulate the operator 
A in terms of operators of Askey- Wilson type, so that we work with a generalized 
hypergeometric operator L. By an Askey- Wilson operator D, we do not only mean 
the operators introduced by Richard Askey and James A. Wilson in [2]. Actually, 
we shall also be working with other kind of operators strongly related with those, 
such as the one introduced by Hahn. Besides, we shall use the companion operator 
S in the sense of Magnus [H]. In particular, the hypergeometric operator we shall 
work with has the form 

L = crD 2 + tSD, 

where a € P2 [x] and r € Pi [x] . The operators S and D are defined by means of a 
two variable quadratic polynomial P which in our case will be symmetric. In our 
context, a translation of the difference operator A/(t) = /(t + 1) — f(t) and the 
classical derivative will be particular cases. These Askey- Wilson type operators act 
over polynomials -also over analytic functions- without any explicit reference to a 
particular sequence Xt satisfying properties (a) and (b). Formulated in these terms, 
it turns out that the operator A becomes more intrinsic and, as we shall see, this 
provides some relevant advantages. 

In order to study the eigenvalue problem posed above, we previously state some 
results for the Askey- Wilson operators. We first give seven canonical forms under 
affine transformations which allow us to provide a very natural classification of the 
Askey- Wilson operators. We also develop the Leibniz and commutation formulas 
for the Askey- Wilson operator and its companion, the corresponding Taylor formula 
and Cauchy's integral formula. Results of this kind are interesting on their own 
right and can be used to obtain new proofs of well-known summation formulas for 
some basic hypergeometric series. We have not included here these applications 
since they will be the subject of a forthcoming publication The use of these 
techniques and mainly our Taylor formula, allow us to obtain simple expressions 
for the eigenfunctions of the hypergeometric operator L in a very general context, 
see Theorem 13.41 In particular, we do not need to invoke any Rodrigues type 
formula or any result from the theory of basic hypergeometric polynomials, such 
as the Watson-Sears transformations used in . Moreover, the identities given in 
Theorem 13 .41 seem to be unknown in such a general setting and, as we shall see, can 
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be easily applied to particular cases. The expressions of the eigenfunctions in terms 
of g-hypergeometric functions arise from this result. On the other hand, Theorem 
13.101 provides an algorithm to obtain Rodrigues type formulas which depend on 
the resolution of certain functional equation. Our methods in Theorem 13.101 are 
closely related to those used in ^1]. However, the final result of Nikiforov, Suslov 
and Uvarov is not so clear since they do not provide the relation with the Askey- 
Wilson operators. Moreover, it is not easy to see the dependence of their Rodrigues 
formulas with respect to the function x t which appears in properties (a) and (b). 
Our Rodrigues formulas are also related to the ones obtained by Askey and Wilson 
in [2] for the Askey- Wilson polynomials. 

The canonical Askey- Wilson operators mentioned above lead us to consider five 
canonical families of eigenvalue problems. For some reasons that will become clear 
throughout the paper, we call these families as follows: continuous, arithmetic, 
quadratic, geometric and trigonometric. The corresponding eigenfunctions include 
as generic cases the Jacobi, Hahn, Wilson, g-Hahn and Askey- Wilson polynomials 
respectively. We shall parameterize each family by means of one or two auxiliary 
polynomials. In the trigonometric case, we get a Laurent polynomial. Using this, 
we rewrite the expressions of Theorem 13.41 in terms of the basic hypergeometric 
functions. Besides, we will study the functional equation of Theorem l3.1(Jl to obtain 
explicit Rodrigues formulas for each canonical case. Our approach looks mainly at 
the regularity of the solutions of the functional equation. This gives rise to new 
Rodrigues type formulas which do not appear in the literature. Quite surprisingly, 
apart from the continuous case, we can always obtain non-ramified solutions. More 
specifically, we get either meromorphic solutions or solutions with one essential 
singularity. In general, there is not uniqueness of solutions, so that we investigate 
a great variety of solutions. Although we shall only give explicit identities for 
some relevant generic cases, it will become clear that our methods can be applied 
in all cases. Moreover, in contrast with ^1], our methods allow us to obtain such 
relations directly. That is, we do not need to take limits. We also point out that our 
Rodrigues formula for the Askey- Wilson polynomials is different from the one given 
by Askey and Wilson in [2J. From our point of view, looking at the regularity of 
the functions appearing in Rodrigues formula, Askey and Wilson obtained ramified 
solutions while we get meromorphic ones. By using the non-uniqueness of the 
associated functional equation, we show how our Rodrigues formula can be used to 
get the orthogonality relations given in In fact, the techniques we develop here 
can be used to obtain discrete and absolutely continuous measures with respect 
to which, the corresponding eigenfunctions become pairwise orthogonal. Besides, 
one of the most relevant features of our techniques is that the role of the theory of 
elliptic functions in the problem becomes quite clear. More concretely, our approach 
divides the analysis of Rodrigues type formulas into two different parts. First, we 
seek for regular Rodrigues type formulas avoiding ramified functions. Formulas of 
this kind are useful for applications in combinatorics. Second, with the aid of the 
theory of elliptic functions, we explore the non-uniqueness of solutions to obtain 
the original Rodrigues type formulas given in 0. 

Throughout this paper, we use the customary notation [H] for ^-shifted factorials 
and basic hypergeometric series. We also use an auxiliary non-zero parameter A 
with q = A 2 . The restriction \q\ < 1 is only assumed in some explicit expressions. 
For instance, the main Theorems 13 . 41 and 13 . 1 01 hold without that restriction. 
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1. P-SEQUENCES 

In this section we introduce the notions of P-sequence and P-function associated 
to a given symmetric polynomial P. As we shall see, these notions become very 
useful in order to give a somehow natural classification of the family of operators of 
Askey- Wilson type. Given a, b, c complex numbers, let us consider the symmetric 
polynomial in two variables P(x, y) — x 2 + y 2 — 2axy — 2b(x + y) + c. For any such 
polynomial we have 

P(x,y) = y 2 -2A(x)y + B(x), 
where A(x) = ax + b and B(x) = x 2 — 2bx + c. Let iZ stand for the set of integer 
and semi-integer numbers. A sequence of complex numbers (xt), with the index t 
running over ^Z, will be called a P-sequence if 

P{xt,y) = (y~ x t+ i)(y - x t _i) for all t S ±Z. 

If (x t ) is a P-sequence and to S ^Z, the translated sequence (yt) with y t = x t -t a is 
also a P-sequence. We shall say that xq is the base point of the P-sequence (xt). 
Given a complex number £, a P-sequence with base point £ can be constructed 
recursively. Moreover, this P-sequence is unique once we have decided which root 
of P(£, y) is X\/2 and which one is x_i/ 2 . In particular, we obtain two P-sequences 
with base point £ which coincide when P(£,y) has a double root. If (xt) is such a 
P-sequence, the other P-sequence with base point £ is (yt) where yt = x-t- 

Remark 1.1. Obviously, the symmetry of P plays an essential role in the iteration. 

As we have seen, P can be regarded as a polynomial in the variable y with 
coefficients depending on x. This allows us to consider the discriminant of P as a 
function of the variable x which, up to a constant factor, has the form 

S(x) = A{x) 2 - B(x) = (a 2 - l)x 2 + 2b(a + l)x + b 2 - c. 
By the definition of P-sequence, the following identities hold for any t € |Z 

A{x t ) = ^(x t+ i+x t _i) 
(!) B {xt) = x t+ ix t _i 

5(x t ) = \{x t+ i-x t _i) 2 . 
Remark 1.2. characterizes P-sequences with (A(xt), B(xt)) or (A(xt), S(xt)). 

Each symmetric polynomial P(x,y) = x 2 +y 2 — 2axy — 2b(x + y) + cis completely 
determined by its coefficients a, b, c. Given a parameterized curve 

7 : K — > C 3 

s i — > (a(s),b(s),c(s)), 

we consider the polynomials Y s (x,y) = x 2 + y 2 — 2a(s)xy — 2b(s)(x + y) + c(s). 
Then, once we have fixed a complex number a^o, we can construct for any s S K a 
P s -sequence (xt(s)) with base point xo such that the whole family of P s -sequences 
depend continuously on the parameter s. Namely, if 6 S stands for the discriminant 
of P s , we can choose the value of \J 8 s {xo) in such a way that we obtain a continuous 
function on s. This allows us to define 

x ±±( s ) = A s(x ) ± y/S s (x ) 

and then we construct the P s -sequence recursively. The continuous dependence 
along curves expressed above shows that any P-sequence with a = ±1 can be 
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approximated by P-sequences with a^±l. This continuity argument will be very 
useful in the sequel and we shall use it in what follows with no further reference. 

Let V be the set of polynomials P(cc, y) = x 2 + y 2 — 2axy — 2b(x + y) + c. The 
set V is an affine subspace of the vector space SP2[x, y] of symmetric polynomials 
in two complex variables having degree < 2. The affine group of the complex plane 
Aff(C) naturally acts on V by the action Ai : Aff(C) x V -> V given by 

(g-P)(x,y)^C 2 P(9~ 1 (x),g~ 1 (y)), 

where g(x) = £x + rj. The parameter a, which will be called the main coefficient 
of P, is invariant under the action of Aff(C) and the discriminant 8 is transformed 
by the following rule 

(2) S*- p (x) = ( 2 8 p (g~\x)). 

The relation above defines another action A2 ■ Aff (C) x P2 [a;] — > P2 [x] on the vector 
space P 2 [x] of complex polynomials in one variable with degree < 2. As we shall 
see throughout the paper, all the notions associated to the elements P of the set V 
are naturally preserved under the action of Aff(C). For instance, if we are given a 
P-sequence (x t ), we have that (g{x t )) is a (g ■ P)-sequence. A P-sequence (x t ) is 
an even function of the variable t if and only if its base point xq is a root of the 
discriminant 6. For this reason we define the even points of P to be the roots of 
its discriminant 6. The number of distinct even points of P will be denoted by evp. 
Now we show that the main coefficient of P and the parameter evp arc relevant 
invariants in the orbits of the action A% . 

Proposition 1.3. Given a pair of symmetric polynomials Pi,P2 € V with main 
coefficients Cij and discriminants 8j for j — 1,2, the following are equivalent: 

(a) The symmetric polynomials Pi and P2 belong to the same Ai-orbit. 

(b) The discriminants 8\ and 82 have the same degree, evp 1 — evp 2 and a\ — 
a 2 - 

Proof. The implication (a b) follows from (J2J and the recalled fact that the main 
coefficient remains invariant under the action of Aff(C). Let us see that (b =>• a). 
If (b) holds, it is clear that Si and 62 belong to the same ^-orbit. In particular, 
there exists f/12 € Aff(C) such that the polynomial P = 512 • Pi has discriminant 82. 
If a ^ — 1 this implies P = P2 and we are done. Finally, for a = — 1 we can always 
choose gi,g 2 & Aff(C) such that 

(ft ;-Pj)(x,y) =x 2 +y 2 + 2xy + C] for J = 1,2 

where Cj = 1 if 8j ^ and Cj = otherwise. Since 81 — if and only if 8 2 — by 
hypothesis, Pi and P2 belong to the same yli-orbit and the proof is completed. □ 

Using Proposition !! ,3l we can write any polynomial P in V, and its corresponding 
P-sequences, in a canonical form. Namely, in the following table we write all 
the ^4i-orbits by considering all the possible combinations of ^-invariants: the 
main coefficient a, the degree dg(8) of the discriminant and the number of distinct 
even points evp. We have chosen a canonical symmetric polynomial P in each 
yli-orbit. Then, we obtain the corresponding P-sequences. Our election for the 
canonical forms has given priority to the simplicity of the resulting P-sequence and 
it is somehow arbitrary. Finally, recalling that Aff(C) can be regarded as the semi- 
direct product 1(C) x GL(C) where 1(C) is the group of translations in the complex 
plane, we include the isotropy subgroup Ip associated to each orbit. In the following 
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table, we use a new parameter A ^ denned by the relation A = a ± \[a?~- ~T. In 
other words, 

A + A" 1 

Obviously, the product of the two values assigned to A is 1. Moreover, we have 

A = 1 ^ a = 1 
A = -1 <^ a = -1 

As we have noticed above, these values of A can be regarded as limiting cases. 





a 


dg(5) 




Canonical form 


P-sequence 


Ip 


T 


¥= ±1 


2 


2 


x 2 + y 2 — 2axy + a 2 — 1 


i(A 2t « + X-^u- 1 ) 


{±1} 


G 


# ±1 


2 


1 


x + y — 2axy 


X ±2t x 


GL(C) 


Q 


1 


1 


1 


(x - y f - i( x + y) + TS 


t 2 + 2t« + u 2 


1 


A 


1 








(* - y) 2 - k 


Xq ± t 


1(C) xi {±1} 


C 


1 


— oo 


oo 


(x - y) 2 




Aff(C) 


O 


-1 








(* + «) 2 - i 


(-l) 2t (x ±t) 


{±1} 


E 


-1 


— oo 


oo 


(x + y) 2 


(-l) 2t xo 


GL(C) 



Table I. Canonical P-sequences. 



The first two canonical forms include infinitely many ^li-orbits, one for each 
value of the main coefficient a. There are only seven inequivalcnt canonical forms. 
The capital letters in the left column are acronyms of the names we have adopted for 
the canonical forms. These names have been chosen attending at the shape of the 
canonical P-sequences: Trigonometric, Geometric, Quadratic, Arithmetic, 
Constant, Oscillating and Even. For some reasons that will become clear later, 
we shall also refer to the Constant form as the Continuous canonical form. 

Remark 1.4. It is clear that the polynomials belonging to the ^4i-orbits associated 
to the trigonometric canonical form are dense in V. Therefore, the trigonometric 
canonical form is the most relevant one since any other canonical form can be 
regarded as a limit of families of polynomials in this particular orbit. 

Remark 1.5. Taking affine transformations from the canonical P-sequences, it is 
not difficult to check that every P-sequence (xt) has one of the following forms: 

• If a ^ ±1, then we have x t = k + k\\ 2t + k 2 X~ 2t where 

x = k + ki + k 2 
b = (1 — a)k 
c = 2bk + Ak 1 k 2 {a 2 -l). 

• If a = 1, then we have x t — ko + kit + k 2 t 2 where 

x = k 
b = \k 2 

C = ig^-2 4^1 kok 2 . 

• If a = —1, then we have x t = k + (— l) 2t (ki + k 2 t) where 

x = k + ki 
b = 2k 
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These relations already appeared in the book ^3] of Nikiforov, Suslov and Uvarov. 

Proposition 1.6. Every V -sequence (x t ) satisfies the following relations, with the 
obvious limits for A = ±1, for any s G 

x t +s + x t -s A 2s + A~ 2s X 2s + \- 2s - 2 , 

x t 1 



2 2 A + A- 1 -2 

A 2s - A- 2s 
A- A- 



X t+S - Xt-s = — TTTf- Ot+i 



Proof. By continuity, it suffices to consider P-sequences with a ^ ±1 . Moreover, 
by a translation in the variable t, we can assume that t — 0. Finally, under these 
new hypothesis, the result follows easily from Remark 1 1.51 □ 

Corollary 1.7. Every P '-sequence (xt) satisfies the following relations: 

(a) x t+ i - x t -! = 2a (x t+ i - x t _i). 

(b) 8a S(xt) = (x t+ i - x t -i){x t+ i - x t _i). 

(c) x t +i - 2x t + x t -i = 4(a + \){A{x t ) - x t ). 

Given a symmetric polynomial P G V, an entire function 9 : C — ► C will be 
called a P-function if the sequence 

...6{z- 1), 8{z - §),«(*), 0(z + \),9{z + 1), . . . 

is a P-sequence for any z e C. To be consistent with our notation, we shall write 
6t for 9(t). P-functions and P-sequences can be treated in the same fashion. For 
instance, an entire function 9 : C — > C is a P-function if and only if the following 
relations hold for any t G C 

A(9 t ) = \{9 t+h + t _i) and S(0 t ) = 9 t+ r0 t _i. 

Remark 1.8. In what follows, when dealing with P-functions, we shall write A* as 
an abbreviation of e tu) for some election uj of the logarithm log A and any complex 
number t. We shall also make this slight abuse of notation for g* where q = A 2 . 

2. Taylor formula for the Askey- Wilson operator 

In this section we develop a Taylor formula for operators of Askey- Wilson type, 
which will we shall apply in Sections [3] and 0] to give explicit expressions for the 
polynomic eigenfunctions of the associated hypergeometric operator. 

2.1. Operators S and D. If ¥[x] stands for the space of complex polynomials in 
one variable, let us denote by Pfc[x] the subspace of polynomials of degree < k. 
Given / G P[x], we define fd(x,y) as the polynomial which coincides with 

f_M - f{y) 
x-y 

whenever x^y. We also introduce the polynomial f s (x, y) — \{f{x) + f{y)). Both 
f s (x,y) and fd{x,y) are symmetric polynomials in two variables. On the other 
hand, given a symmetric polynomial P(x, y) — x 2 + y 2 — 2axy — 2b(x + y) + c, we 
shall denote by u(x) and v(x) the roots of P(x,y) in the variable y. That is, we 
have P(x, y) = (y — u(x))(y — v(x)) for all ieC. Then we define the operators S 
and D as follows 

S/(x) = f s (u(x),v(x)) and D/(x) = f d (u(x),v(x)). 
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As it is recalled by Magnus the divided difference operator D is essentially 
the Askey- Wilson operator while S will be called the companion operator to D. 
Any symmetric polynomial in two variables provides a polynomial in the variable 
x when it is evaluated at (it(x), u(x)). In particular, we have that S/ £ P[x] and 
D/ £ P[x] for any / £ P[x]. That is, S,D e End(P[x]) are linear mappings in P[x]. 

Remark 2.1. S and D are naturally transformed under the action of Aff(C). That 
is, if S p and D p stand for the operators associated to the symmetric polynomial P 
and we take some g — (n, () £ 1(C) x GL(C) = Aff(C), then we have 

g S p g- 1 = S g P and g D p g' 1 = ( D s P 

where we identify g £ Aff(C) with the induced mapping (g ■ f)(x) — f(g (x)). 

Proposition 2.2. We have S(P fc [x]) C P fe [x] and D(P fe [x]) C P k -i[x] for k > 1. 
Moreover, if f(x) = x k we obtain the following expressions 

sm = • ... 

A fc - \- k 
X-X- 

where the dots stand for terms of lower degree. If f £ ¥q[x], Sf = f and Df = 0. 

Proof. The action of S and D on Po[x] is trivial. Now, let /(x) = x k for some k > 1. 
Then, since the roots of P(x, y) are given by A(x) ± y/ 6(x), we have 



Df( , = (A(x) + ^6(x))-(A(x)-^5(x)y 
whenever 5(x) ^ 0. The binomial theorem gives 

B fW = ^(^Mx^i^stxjy-^i-i-iy) 

j=0 



2vi* . ■ . 

3=0 

= - f' 1 [(a + ^~l) k -{a- Va^l) k ] +... 
2\/a z — 1 

where the dots stand for terms of lower degree. This new expression is now valid 
for any x £ C. In particular, we have obtained 

\k \— k 

D/(x) = ^-A_x fc - 1 +... 

Similar arguments provide the expression given above for S/(x). □ 

Remark 2.3. Both operators S and D depend polynomially on P. In particular, 
Proposition 12.21 and the forthcoming expressions hold for A = ±1 by taking the 
obvious limits. This follows from a simple continuity argument. 

Proposition 2.4. S and D satisfy the following Leibniz rules 

S(fg)(x) = Sf(x)Sg(x)+S(x)T>f(x)T>g(x) 

D(/ 5 )(x) = D/(x)S 5 (x) + S/(x)D 5 (x). 



CLASSICAL HYPERGEOMETRIC POLYNOMIALS 



9 



Proof. Both expressions are straightforward and we leave them to the reader. □ 

If (xt) is a P-sequence, then we can use our definitions of the operators S and D 
to obtain the following relations 

at, n ,/ x f(x t+ i 2 ) + f(x t _i) 
Sf(x t ) = f s (x t+ i,x t _i) = 2 — 2 — 

Df(x t ) = f d {x t+h ,x t _.) = Xt+h - Xt _ k2 ■ 

Also, it is not difficult to check the validity of the following identities 

S/(zt) = f(x t+ i) - Xt+h ~ Xf -* D/(g t ) 
Sf(x t ) = f(x t _ i )+ +5 2 S- D/fe)- 

Proposition 2.5. S and D satisfy the following relations 

DSf(x) = {a + l)(A(x)-x)D 2 f{x)+aSDf(x) 
S 2 f(x) = ad(x)D 2 f(x) + {a + l){A(x)-x)ST)f(x) + f(x). 

Proof. Given a complex number x, let (&t) be a P-sequence with base point x. By 
continuity we can assume that 

x t+ i^x t _i for t = 0,±f. 
Under these assumptions, we can use Corollary 1 1 . 71 and to obtain 

DS/(x) = i + f(x) + f(x^ y 

xi — x_i V 2 2 

2 2 

Xi — X l \ Xi — X X — X_i / 



(a+i)^)-,)^/^) = ^-te+t-ufw-fw . fw-n*-!) ) 

4(xi — x_i ) \ x\ — x x — X-i J 

2 2 

(a + l)(A(x)-x)SDf(x) = + + /(*) - /(*-ih 

8 V xi — x x — X—i / 

8 V si — x x — x_i / 



S 2 /(x) 



8 V xi — x x — x_i 

/(xr) + /(x) , f(x)+f(x-i) 



4 4 

At this point it is easy to check the stated relations. This completes the proof. □ 

Remark 2.6. Once we have introduced the operators S and D, we are in position 
to justify some terminology introduced above. It follows from Table I that the 
P-sequences of the canonical form C are constant. In particular, S is the identity 
operator while D is given by the classical derivative 

ax 

This is why we have decided to call this canonical form the continuous form. On 
the other hand, the P-sequences of the arithmetic canonical form are given by 



10 



MARCO AND PAR.CET 



xt = xo ± t. In particular, S can be regarded as an arithmetic mean while D can 
be rewritten in terms of the classical operator Af(x) = f(x + 1) — f{x) as follows 

In summary, these canonical forms can be regarded as the 'classical' forms. 

Remark 2.7. By Proposition ^. 51 we can characterize those polynomials P for which 
we have SD = DS. Namely, the operators S and D commute if and only if the 
coefficients of P satisfy a = 1 and b — 0. In other words, this happens only for the 
classical forms described in Remark [ 



2.2. Taylor formula. In order to construct the polynomials which will appear 
in the Taylor formula for the Askey- Wilson operator, we shall need the following 
lemma, which is a simple consequence of Proposition 

Lemma 2.8. Given a P -sequence (x t ), let us consider the sequence yt = x 3 t for 
some integer number j £ Z. Then we have 

- = a i y t + b j and (y t+ i - y t _i) = 4 { x _ X -i ) S W) 

where the coefficients a,j and bj are given by 

X> + \- j , , X 3 ' + \- j - 2 , 

a., = and bj = ; b. 

3 2 3 A + A~ 1 -2 

It turns out by Lemma r2.8l that (yt) is a Pj-sequence, where Pj is the following 
symmetric polynomial 

/A J — A~ J \ 2 f\ ] — A~ J \ 2 

p j(x,y) = (y~ ajx - bj) 2 - { x _ x -! J S ( x ) with S j( x ) = \ x _ x -i ) s ( x )- 

The symmetric polynomials Pj are naturally transformed under the action of Aff(C). 
Namely, we have 

(9-P)j(x,y) = (g-Pj)(x,y). 
Given a positive integer k > 1, we define the polynomials $k(x,y) as follows 

fe k 

®2k(x,y) = Y[P 2 j-i(x,y) and <$> 2 k+i(x,y) = (y - x) Psj(x,y). 

3=1 3=1 

We set $o(x,y) = 1 and <&i(x,y) = y — x. In particular, $k(%,y) = ( — l) fe $fc(y> 
for all k > 0. Moreover, given two P-sequences (xt) and (yt), it follows by Lemma 
EH that 

fe-i fe-i 

(5) $fcOo, Vo) = Y\_(Vo - x j-h^l) = Y\_(yj-h^± - x o) for k > 1. 

3=0 " j'=0 

Remark 2.9. The polynomials depend continuously on the coefficients a, 6, a In 
particular, once again we can take limits A — ► ±1 in the forthcoming expressions. 

Lemma 2.10. The following relations hold for any k > 1 

X — k X k 

D#jk(-,p) = A _ A-1 (-,!/) 

D$ fe (av) = x _ x -i *fc-i(ss,0- 
Moreover, for k — we have the relations D<I>o(-, y) = D$o(£, •) = 0. 
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Proof. Given a complex number x, we consider a P-sequence with base point x. 
Moreover, by continuity we can assume that X1/2 7^ X-i/2- In that case, it is not 
difficult to check that 

$fc0i/2,y) - ®k O-1/2, y) a;-fe/2 - ^fc/2 _ , . 
= 9 k -i{x ,y). 

^l/2~ x -X/2 a; l/2 — ^-1/2 

The first relation is then obtained by Proposition The second relation follows 
by the first relation and the identity $fc(x, y) = (— l) fc <&/.(y, x). Finally, when k = 
the stated relations are obvious. This completes the proof. □ 



Lemma 2.11. The linear operator dk ■ ^d[x] — * P<2_ fc[x] defined by 

A -A' 1 

\k-j _ \j 



d «fW = ( H TI^-A-t) D "/0) Aw |A| * 1 



/ias a continuous extension for any A which depends polynomially on P. 

Proof. From Lemma 12.101 the following relation holds whenever r > k and the 
parameter A satisfies |A| 7^ 1 

fc— 1 

(6) wx, .) = ( n A x l ; p I') •). 

Moreover, 9fe$ r (x, •) vanishes for r < k. By a well-known result of g-combinatorics, 
see for instance |1U| . the rational function which appears in JIjJl is a symmetric 
polynomial in A and A -1 . In particular, the function <9/c$ r (x, •) can be written 
as a polynomial in a = |(A + A _1 ),6 and c. Since for any fixed x, the family 
<&o(x, •); $iO&j ')>••■ ^d{x, •) is a basis of Pdfx], the proof is completed. □ 

Remark 2.12. The operator 9ft defined in Lemma r2.11l can be regarded as the analog 
of the divided power of the k-th derivative d^- k \ which appears in finite field theory. 

Remark 2.13. Taking q — A 2 , the rational function appearing in the proof of Lemma 
12.111 can be rewritten in terms of a g-binomial coefficient 



fc-i 



A r 3 - X> r k{r _ k , 

llA fe -^-A^- fe 

3=0 



where 



1 {q;q)k(q;q) r -k 



The following result constitutes a Taylor formula for the Askey- Wilson operator. 
As we shall see immediately, it generalizes the continuous and discrete classical 
expressions of this formula. 

Theorem 2.14. If (xt) is a ¥ -sequence and f £ P r [x], then we have 

r k— 1 

f(y) = J2 d kf( x k/2) Y[(y-xj). 

k=0 j=0 

Proof. There exists a family of complex coefficients Ao, Ai, . . . A r such that 

r fc— 1 

f(y) = J2 x « n^-^)- 

fe=0 3=0 
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fc-1 

On the other hand, since (y — Xj) — ^fc(^(fc-i)/2j2/)i formula © gives 



k=j 1=0 

By identity 0, if k > j then <5>fc_j(x(fc_i)/2, y) vanishes at 

y — x i+ i for i = 0,l,...k — j — l. 

In particular, evaluating our expression for djf(y) at the point y — Xj/ 2 , we obtain 
the identity A& = dk.f{xk/2) as we wanted. This completes the proof. □ 

Remark 2.15. The classical continuous and discrete Taylor formulas are particular 
cases of Theorem 12.141 Namely, if P(x,y) — x 2 + y 2 — 2xy the P-sequences are 
constant and the operator D coincides with the classical derivative, as it was recalled 
in Remark 12.61 Therefore, Theorem 12 . 141 takes the classical form 

k=0 

On the other hand, if P{x, y) = (x — y) 2 — j, we know that the P-sequences have 
the form xo ± t. Thus, using the classical notation [t]k = t(t — 1) • • • (t — k + 1) and 
Af(x) = f(x + 1) — f(x), we obviously have 

1 

fc! 

In particular, Theorem 12 . 1 41 provides the discrete Taylor formula 

k=0 ^ 

2.3. Some remarks for analytic functions. Given an open subset fl of the 
complex plane, we denote the space of analytic functions in Q by 7i(f2). Also, 
7 ~ (modfi) means that 7 is a cycle in homologous to zero with respect to SI. 
Finally, given z G Q, Ind(7, z) denotes the index of z with respect to 7. 

Lemma 2.16. Given an open subset f2 of the complex plane, let us denote by A.q 
the diagonal of S! x fl. Then, for any f E 7i(Sl), the function 

f d (u,v)= IM^IM f or ( u ,v)e(nxQ)\A Q 
u — v 

can be continuously extended to an analytic function : f2 x 11 — > C. Moreover, if 
7 ~ (modi!) and Ind(7, u) — Ind(7, v) — 1, then we have 

fd(u,v) = — / 7 T7 zdy. 

2tti J 1 [y -u)(y- v) 

Proof. It is a simple consequence of Cauchy's integral formula. □ 



$ k (x,y) = [y-x+^] k and d k f(x) = - A k f(x - § ) 



Let us consider the set fi^ = {x £ C : A(x) ± 8{x) £ fl}. By Lemma 12.161 
we can define the operator D : 7i(Sl) — * W(fiW) as follows 

1 r m 



(7) Df(x) = f d (A(x) + y/8@j, A(x) - y/8(xj) = —J dy, 
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with 7 ~ (mod O) and Ind(7, A(x) ± \J 8{x)) = 1. It is clear that extends the 
original definition of the Askey- Wilson operator D. Moreover, we can also extend 
the definition of the companion operator S. Namely, if / £ H(f2), xo £ and 
S(xq) ^ 0, the function 



S/(a;) = f s {A(x) + y/S^,A(x)-y/6(^) = 



f{A(x) + ^6{x)) + f(A(x) - y/8ffi) 



2 

is obviously analytic for a; in a neighborhood of xq. Besides, S/ is continuous at xq 
when S(xq) = 0. Thus, the roots of 5 are removable singularities and Sf becomes 
analytic in f^ 1 ). If we define recursively the sets 

0( fc+1 ) ={ieC: A(x) ± y/S(xj £ n {k) } with ft (0) = f2, 

we can consider the iterated operators S fc ,D fc : 7i(0) — * 7i(Sl( fc )). An open subset 
of the complex plane will be called P-invariant if, for any P-sequence with base 
point xo £ f2, we have x±i/2 £ ^- In other words, if fi C Q^' for any positive 
integer k > 1. If O is P-invariant, we deduce that S fc / £ W(f2) and D k f £ H{Q) 
for any function / analytic in Q and any positive integer k. 

Remark 2.17. Similarly, if .M(f2) stands for the space of meromorphic functions in 
f2, we can extend the previous operators so that 

S fc ,D fe : M(n) -> M(n ik) ). 

Furthermore, if is P-invariant, then we have operators S fc ,D fc : .M(f2) — > A4(Sl). 

Remark 2.18. Obviously, ProDOsitions l2.4l and l2.5l remain valid in this new context. 

Proposition 2.19. Let fi be an open set of the complex plane and 7 ~ (mod Q). 
Then, given f £ Tl{VL), an integer k > and x £ fl^ k \ we have 

if x is the base point of a P-sequence (x t ) with Ind(7, Xj_k) — 1 for j — 0,1, ... k. 

Proof. The cases k = and k — 1 are Cauchy's integral formula and (0 respectively. 
The general case follows by induction from J2J and the relation 

A fe - \- k 



D 



(a) = A _ A _i $k+i(x,y) 1 



-*k(-,tf)- 

which can be checked by the reader. This completes the proof. □ 

Corollary 2.20. Let be an open subset of the complex plane and 7 ~ (mod fi). 
Then, given f £ H(Q,), the expression 

dk ^ = Y--l^r-^ 

2m J 1 9 k+ i(x,y) 
defines a linear operator d k ■ H(£l) — > 7i(il^ k '). Moreover, we have 

Remark 2.21. It turns out that Lemma l2.11l is now a consequence of Corollarv l2.2UI 
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Remark 2.22. When x is the base point of a P-sequence (xt) with x^_k pairwise 
distinct for j — 0, 1, . . . k, we have by Corollary 12. 201 



W = En r- 

j=0 lltKi^Xfc^j 



As we explain in Remark 13.111 this is useful to display our Rodrigues formulas 
given by Theorem 13.101 We recall that Nikiforov, Suslov and Uvarov give similar 
expressions for less intrinsic operators in 14, 3.2.3]. 

Remark 2.23. Under the hypothesis of Remark 12. 2 21 Theorem 12. 141 gives 

k=0 j=Q lll 7=J v ■> % > j=0 

which is Newton's divided difference formula for the interpolation polynomial. In 
fact, as the referee of this paper has pointed out, the formula given in Remark l2.22l 
can be proved by induction. In particular, it is possible to obtain an alternative 
proof of Theorem 12 . 1 41 from the theory of interpolation polynomials. We prefer our 
proof since it shows how the polynomials are related to D, see Lemma \'2. 101 

Remark 2.24. The operators dk '■ Tt(fl) — » Ti(n^) provide the natural framework 
to develop a generalization of Theorem 12.141 Namely, in a separate work ^3] > we 
study the convergence of the Taylor series associated to the operators dk of a given 
function / S 7Y(0) and its relation with the basic hypergeometric functions. 

3. Taylor coefficients and Rodrigues formula 

In this section we introduce the corresponding hypergeometric operator. Then 
we compute the Taylor coefficients of its eigenfunctions. Also, we provide a general 
procedure to obtain Rodrigues type formulas. At the end of this section, we shall 
study the discrete orthogonal relations which arise from our techniques. 

3.1. The hypergeometric operator. Given a symmetric polynomial P in the set 
V, we consider the associated Askey- Wilson operator D and its companion S. Let 
a E P-2 [x] and t € Pi [x] given by 

cr{x) = CC2X 2 + ol\X + ao and t(x) — [3\X + /3 . 

The aim of this section is to study the polynomic eigenfunctions of the following 
operator 

L = o-D 2 + tSD. 

L will be called the hypergeometric operator associated to P, a and r. In 
Section 0] we shall study in detail some relevant particular cases. It follows by 
Proposition 12 . 21 that L(Pfe[x]) C Pk[x] for k > 1. Moreover, if f(x) ~ x we obtain 

Lf(x) = -fi k x k + ... 

where the dots stand for terms of lower degree and 

A fe - A k \( A- A 1 „ \ k -i (' A- A 1 \ j_ fc 



a 2 + A A - a 2 ^ /3i A 



(A -A- 1 ) 2 

As an endomorphism of [x] , the hypergeometric operator L depends continuously 
on P, er and r. Therefore, taking suitable limits for A = ±1, it turns out that ^ is 
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a continuous function on A. The following formula, for which we should take limits 
in A when A = ±1, will be useful in Section 0] 



A — A n \ , ,-j_t._i / A — A 

(A- A- 1 ) 2 

Now we point out some relevant remarks on the operator L: 

• If L p stands for the hypergeometric operator associated to the symmetric 
polynomial P and we consider some g = (77, Q £ 'X(C) x GL(C) = Aff(C), 
it follows easily from Remark 12. II that 

g LV 1 = CV(D 9 ' P ) 2 + Cr 9 S 9 P D9 p 

where (g ■ f)(x) = f(g~ 1 (x)), a 9 {x) = a(g^ 1 (x)) and t 9 (x) = r(g^ 1 (x)). 

• When (xt) is a P-sequence with x t _i 7^ x t+± f° r i = 0, ±i, we can write 

L<f{XQ) = 



Xl/2 - X-!/ 2 



Xi — Xo Xq — X-i 



T~(x )\f(Xi)-f(x ) , f(x ) - 



Xi — Xq Xo ~ X-i 

In fact, there exist a limiting version of this equality for any P-sequence (xt). 
Moreover, since x\ — x_\ = 2a(xi/ 2 — £-1/2) by Corollary 1 1.71 L/(xo) can 
be written in terms of Xt with t £ 7L. In other words, since the subsequence 
of (xt) which arise when t runs over Z is a P2-sequence with the terminology 
introduced after Lemma 12.81 we can say that L only depends on P2. 
If the symmetric polynomial P is of type O (resp. E), it can be checked 
that P2 is of type A (resp. C). Besides, if P is of type A (resp. C), we 
also have that P2 is of type A (resp. C). Therefore, taking into account 
the fact that L only depends on P 2 , we do not loose generality if we restrict 
our study to the canonical forms T, G, Q, A and C. This point of view 
leads to a different and somehow more natural classification than the one 
given in ^] 3.4]. Namely, while our classification is purely complex, it can 
be said that the classification given in lives over the real field. As we 
point out below, another remarkable advantage of our formulation is that 
affine transformations become trivial. 

If A = ±i, we have D 2 = (A + A -1 )^ = 0. In particular, L = rSD and so 
we get a degenerate case. This phenomenon can be avoided by considering 
the more general operator 

L' = ad 2 + rSD. 

It is not difficult to adapt the methods we shall develop to this case. When 
A = ±i in P, it can be checked that P2 is of type E. Therefore, the even 
canonical form can be regarded from this alternative point of view. 
The polynomial P2 can not be of type O. However, we can construct a 
limiting operator which somehow corresponds to certain P2 of type O. 
Namely, it suffices to take a sequence P(n) in V such that the sequence 
P(n)2 converges in V to a symmetric polynomial of type O. It turns out 
that, in this situation, the sequence P(n) does not converge when n — > 
00. However, it can be checked that the corresponding hypergeometric 
operators L(n) do converge. Therefore, the oscillating canonical form could 
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be analyzed by means of the operator 

L = lim L(n). 

n — >oo 

Although it can be studied what happens with our methods after taking 
this limit, we shall not cover this problem here. 

3.2. Taylor coefficients and Rodrigues formula. The first step in our process 
is the following commutation relation for the iterated hypergeometric operators. 
Expressions of the same kind can also be found in |2j . 

Lemma 3.1. Taking a® = a , tq = r and given j > 0, let us define the polynomials 
a j G P2 [x] and tj € Pi [x] inductively as follows 

aj + i(x) = So-j(x) + (a + l)(A(x) — x)STj(x) + a5(x)T)Tj(x) 

Tj + i(x) = ~D<jj(x) + aSrj(x) + (a + l)(A(x) — x)~Dtj(x) 

Then, the iterated hypergeometric operator hj — cr^D 2 + TjST) satisfies 

D ( L j - Mi) = ( L j+i - Mj+iP- 

Proof. Applying the Leibniz rules and the commuting relations developed for the 
operators S and D in Propositions 12 . 41 and 12.51 we easily obtain 

DL,/(x) - D^-DVX^+D^SD/)^) 

= S(7j(.t)D 3 /(>) + Ba 3 (x)SI) 2 f(x) + Sr i (a:)DSD/(x) + BT j {x)S 2 Bf(x) 

= S<T j {x)D 3 f(x)+Dcr ] {x)ST) 2 f{x) 

+ Sr i (x)((a + l)(A(x) - x)D 3 f(x) + aSD 2 /(x)) 

+ DT 3 (x)(aS(x)D 3 f{x) + (a + l){A(x) - x)SD 2 f(x) + D/(x)) 

= a J+1 (x)B 3 f(x) + r J+1 (x)SB 2 f(x) + D Tj (x)Df(x). 

In other words, T)Ljf(x) = (j^j+i + Dtj)D/(x) where Dr., is a constant Vj since 
Tj 6 Pi [a;]. On the other hand, taking f(x) = x k , we define fij^ by the relation 

Ljf(x) = -n 3: kx k H 

where the dots stand for terms of lower degree. Then, Proposition 12 . 21 gives 

y ^k k v / ^k ^ k 

-(y^-rjMM- = (- — )•: >>,-:.>. 1 • <V- 

By continuity, we have Atj+i.fe— 1 = (ij,k + v j or equivalently 

i-i 

f-j.k = f-j + k + I'.- 

i=0 

In particular, we obtain that p,j + vq + . . . + Vj—i = fij,o = 0. Therefore, we deduce 
the identity Vj = fj,j — /ij+i as we wanted. This completes the proof. □ 

Given a symmetric polynomial P(x, y) - - x 2 + y? — 2axy — 2b(x + y) + c, let us 
take a P-function Qt- Then we define the auxiliary functions 

, e t+ i - e t _ x 

af(t)=a j (9 t )±^- 2 *- Tj (9 t ) for j > 0. 

We shall write o ;± to denote the functions respectively. 
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Lemma 3.2. The auxiliary junctions satisfy the following relations 

af{t)=a+(t+t) and aj (t) = c?-(t - §). 
Proof. By definition we have 

= S^(0 t ) + (a + l)(A(8 t ) - t )STj(0t) + a5{e t )VT 3 {6 t ) 

± \{9 t+h - t _*)(D^(0 t ) + oSr^flt) + (a + l)(A(0 t ) - e t )D Tj (9 t )). 

Then, by Corollary 1 1 . 71 and we get 

~± , + s + t+1 -2O t + t - l f . 
°f+i(*) = ~ + § (r,(^ + i)+r,(^_i)) 

+ ^^^ (rj^-r^)) ± ^±i^i( Ti (^) +r,(^_i)) 

± g {T j (6 t+ t)-T j (6 t _i))± 2 — '-. 

Simplifying the expression above, we obtain &f + i(t) = &f(t ± ^). Therefore, the 
stated relations follow easily by induction and so the proof is completed. □ 

Now, given a positive integer d, we assume that jjij ^ fik for < j < k < d. 
In particular, this condition implies that L is a diagonalizable operator in P<j[a;]. 
Therefore, for any < k < d, there exists an eigenfunction fk € Pfc[x] with degree 
k satisfying 

Lfk(x) + fJ,kfk(x) = 0. 
Moreover, if we define fjf ^ — D J /fc, then Lemma fo. II and induction gives 

(8) ^^ + (^-^=0. 

Lemma 3.3. The following relation holds for any j > 

zj(t)fj? +2 \6 t ) + T 3 {6 t )fl? +1 \e t+h ) + ( W - = o. 

Proof. Applying the definition of aj in terms of a 3 and tj and relation P|. we can 
rewrite the left hand side as 

lhs = vjW&fjpVt) + Tj (e t )SDfi j \e t ) + - ^)fi f) (o t ). 

In other words, the left hand side coincides with Ljfu\0t) + (fJ-k — Mi)/fc \®t)i 
which vanishes by relation JHJ. Therefore, the proof is completed. □ 

If we take the P-function 9 t in such a way that a~ (0) = 0, then we can give 
explicit expressions for the eigenfunctions /o, /i, . . . fa in terms of the Taylor formula 
developed in Section 

Theorem 3.4. Let us assume that <x - (0) = ; then the eigenfunctions fo, fi, • ■ ■ fd 

of L satisfy the following identities 



Aw-*A£(n' A1 t:*:r^ )n(--») 
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where dkfk stands for the coefficient ofx k in fk(x). Moreover, assuming fk(0o) 7^ 0, 
we also have the following identities 

aw - aw t ( n n^^i s^y) n' (i - *>■ 

j=0 i=0 v ' ' i=0 

Proof. Since 5~(0) = 0, we obtain from Lemmas 13 . 21 and 13 . 31 the identity 

Jfe (^/2j - —/fe ^(j + l)/2j- 

If we write this recurrence in terms of the divided derivative operator, we get 

8f(f) s_ A^-A^ r 3 (g 3/2 ) 

- A _ A -i 77—^ ^ +1 / fc (0 (j+1)/2 ). 

Then, the first relation follows by iterating this recurrence and applying Theorem 
12.141 The second relation is an obvious consequence of the first relation and the 
identity 

' 1 A 1+i - A- 1 -* 7-i(0 i/2 ) 



v- 1 

which follows from the expression above for j = 0. This completes the proof. □ 



A - A 1 Hi- Hk 



Remark 3.5. Although using a less clear notation, Nikiforov, Suslov and Uvarov 
gave in ^] 3.1.27] an equivalent identity to Lemma 13.21 Moreover, Lemma 13.31 
is quite close to 3.1.24]. However, since there is no Taylor formula for the 
Askey- Wilson operators in , it seems that the authors did not see the relevance 
of these results to obtain explicit polynomic expressions for the eigenfunctions. 

Remark 3.6. Replacing t by the P-function 6' t — 0±ft-t o )> we can use Theorem 13.41 
to obtain explicit expressions of the eigenfunctions fk for each to satisfying either 
<? + (^o) = or (J~(<o) = 0. Such points to are those for which t „ is a root of the 
polynomial Q(x) = u(x) 2 — 5(x)t(x) 2 . Since the degree of Q is less or equal than 
4, we obtain in the generic cases four explicit expressions. We shall see in Section^] 
that there are only two cases, with distinct eigenvalues, for which Q has no roots. 

Remark 3.7. The assumption fij ^ jik for < j < k < d on the eigenvalues is too 
restrictive in Theorem 13.41 Namely, assuming that <? - (0) = 0, it can be proved 
that L is diagonalizable in ¥d[x] if and only if 



' 1 A!+* - A- 



n 



A - A- 1 



i(d i/2 ) = 



whenever fij = fik for some < j < k < d. In this case, the first identity of 
Theorem 13.41 holds summing from jf. = 1 + max{j : fj, } - = - /ifc}. In any case, for the 
sake of clarity, we shall assume that the eigenvalues are pairwise distinct. 

The first relation in Theorem 13.41 is the appropriate one when we normalize the 
eigenfunctions fk so that they become monic polynomials. The second relation in 
Theorem 13.41 will be very useful in Section since it will allow us to express the 
eigenfunctions fk in terms of basic hypergeometric functions. The next result can 
be regarded as an algorithm to obtain Rodrigues type formulas in this setting. We 
recall that the P-invariant domains were already introduced in Paragraph s. 31 
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Lemma 3.8. Let p and pi be meromorphic functions in for some P -invariant 
open set f2. Then, the following are equivalent: 

(a) For any function f 6 M(fl), we have 

p(x)Lf(x)=D( Pl Df)(x). 

(b) The functions p and p± satisfy the relations 

p{x)a{x) = Sp\{x) and p(x)r(x) = Dpi(x). 

(c) The functions p and p\ satisfy the relations 

(cl) D(pa)(x) = (a + l)(A(x) - x)D(pr)(x) + aS{pr)(x). 

(c2) pi(x) = S(pa)(x) - a5{x)D{ P T){x) - {a + l)(A(x) - x)S(pr){x). 

Moreover, if the P-function Ot is not 1-periodic, they are also equivalent to 

(d) The functions p and p\ satisfy the relations 

(dl) p{6 t _ k )a+{t - |) = p{6 t+k )d-{t + I). 

(d2) Px{6 t ) = \{p{6 t _0a+{t-\)+p{6 t+h )a-{t+\)). 

Besides, the functional equation (cl) (resp. (c2)) is equivalent to (dl) (resp. (d2)). 
Finally, if (a) holds then pi satisfies (cl) and (dl) for the iterated operator Li. 

Proof. Properties (a) and (b) are equivalent since, by Proposition ^. 41 we have 

p(x)hf(x) = p(x)a(x)B 2 f(x)+p(x)T(x)SBf(x) 

D( Pl Df)(x) = (Sp 1 )(x)B 2 f(x) + (Dp 1 )(x)SBf(x). 

To see (b c), we apply Proposition 12. 51 

D(pa)(x) = {a + l){A{x)-x)D 2 p 1 +aSDp 1 (x) 

S(pa)(x) = a6(x)T>( P T)(x) + (a + l)(A(x)-x)S(pT)(x)+p 1 (x). 

Let us denote by L(cl) (resp. R(cl)) the left (resp. right) hand side of (cl). For the 
equivalence between (cl) and (dl), we use Corollary 1 1.71 and Namely, arguing 
as in Proposition 12. 51 we easily get 

cr-(t+±)p(O t+ i)-Z+(t-±)p(O t _i) 

The equivalence between (c2) and (d2) follows by similar arguments. Now, let us 
assume that property (d) holds. Then we will have 

Pl (e t+i ) = p(e t )v+(t) = P {e t )(a{e t ) + ° t+i ~ ±3 T (9 t j) 

Pl (9 t _ h ) = P (6 t )a-{t) = p{e t )(a{9 t ) ^ ~ ±* r(9 t j) . 

In particular, Spx(0t) — p{9t)v{9t) and ~Dp\{6 t ) = p(9 t )T(9 t ). That is, we have seen 
(d =>• b). It remains to see that (a) implies that p\ satisfies (dl) for Li. Since (a) 
implies (d), we have by Lemma 13.21 

Pl (9 t _i )a+(t - i) = P (9 t )a-(t)a+(t) = P x{6 t+ ,)a^{t + \). 
But this relation is precisely (dl) for Li. Therefore, the proof is concluded. □ 
Remark 3.9. (dl => cl) is the only implication that fails for 1-periodic P-functions. 



(L(cl)-R(cl))(0 
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Theorem 3.10. Let p 6 A4(fl) for some P -invariant open set O and assume that 
p satisfies the functional equation 

p(9 t )a+(t) = p(e t+1 )a-(t+l). 

Then there exists a family of functions pj £ M.(Q), with j > and such that po = p, 
determined by any of the identities 

Pj+1 (6 t ) = Pj (e t _ i )v+(t + &p), 
p j+1 (O t ) = Pj (e t+h )a-{t-^). 

Moreover, this family of functions provide the following Rodrigues formula for the 
eigenfunctions /o, /i, • • • /d of L 

( n> - nk)) P {x)Mx) = d k f k (fl A ' A J _~ A AJ i h )v k Pk (x). 

j=o j=o 

Proof. The recurrence relations given above for pj follow from Lcmmas l3.2l and l3.8l 
Besides, by Lemma we have pj(x)Ljf(x) = D(pj+iDf)(x) for any / e M.(Ci). 
If we combine this relation with (JSJ, we get 

(p, - p k )p J (x)fi j \x) = D(p j+1 fi j+ V)(x). 
Then we differentiate successively this expression to obtain 

( n>* - ^)) P {x)h(x) = B k ( Pk fi k) ) = d k f k (ji xk ~ j ~*~ k )v k Pk (x) 

where we have applied Proposition 12. 21 This completes the proof. □ 

Remark 3.11. Let us recall that, assuming that x is the base point of a P-sequence 
with Xj_k pairwise distinct for j = 0, 1, ... k and with the aid of Remark 12. 221 we 
can display Rodrigues formula as follows 



D fc p fe (x)= (n 



' 1 \ k - j ~ X j - k \ J-y Pk(Xj 



j=0 j=0 '■"^3 

Remark 3.12. The recurrence relations for the functions pj given in Theorem 13. 101 
have the corresponding analogs in |141 3.2.11 and Pg 64]. However, the authors of 
[T4] work with the functions P^ SV defined by 

pf U (*)=P# s+ |)- 

Hence, these functions can not be easily regarded as functions of the variable 8 S . 
One of our main contributions is that we show how to find meromorphic solutions p 
and pj (with the exception of the form C, see Section 4) in the variable 9 S . To that 
aim, it is essential to regard these recurrence relations in its most intrinsic form as 
in (c2) of Lemma f3. 81 This result has no analog in [14] . 

Remark 3.13. Another remarkable advantage of this alternative formulation is that 
affine transformations become trivial. A good example to illustrate this is Fischer's 
paper .7 j, where the author gets a Rodrigues formula for the q-Racah polynomials. 
Fischer combines results from and |14j with some independent computations. 
However, from our point of view, these kind of results become trivial consequences. 
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Namely, it is very well-known that the q-Racah polynomials R& are related to the 
Askey- Wilson polynomials P& via 

Rk(x; a, f3, 7, 5) = V k (x; a, b, c, d) 

where x = 2^/q r y5 x and 



= VqjS b = ay/ q/jS 

c = f3y/Sqpf d = y/jq/6. 

Therefore, if we are given a Rodrigues formula for the Askey- Wilson polynomials 
(see Section 0} and we want a Rodrigues formula for the g-Racah polynomials, it 
suffices to take /?fc(X) = Pk{x) and substitute the operator D, associated to the Pfc's, 
by the operator D associated to the Rfc's. 

Remark 3.14. The version of Theorem 13 . 1 Ul for the continuous canonical form, for 
which every P-sequence is constant and hence 1-periodic, is very well-known. In 
this case we have <Jj(x) = cr(x) and Tj(x) = t(x) +ja'(x). Moreover, the functional 
equations given in Theorem 13.101 must be replaced by 

f (x) — p(x)t(x) and Pj(x) — p(x)aHx). 
dx 

This also follows from Lemmas 13 . II and 13 . 81 But a direct proof is quite simpler. 

3.3. Discrete orthogonality relations. In this paragraph we shall assume the 
existence of mcromorphic functions p and pj satisfying the hypothesis of Theorem 
13. 101 and we show how these functions can be used to obtain discrete orthogonality 
relations. Given /, g G F[x], wc define h(u,v) to be the symmetric polynomial 
which coincides with 

f{u)g(v) - f{v)g(u) 
v — u 

whenever u ^ v. Then, we define the Wronskian W : V[x] x F[x] — > F[x] to be the 
bilinear form 



W(f,g)(x) = h(A(x) + y/5(x),Mx) - y/5(x))- 
Recall that W(f,g) £ P[x] for any f,g£ ¥[x] since h(u,v) is symmetric in (u,v). 

Lemma 3.15. We have p(x)(fLg — gLf)(x) = D(p{W(f, g))(x) for all x G Q. 

Proof. Given a point x £ J7, we know by Lemma 13.81 that 

p{x){fhg - ghf){x) = (fB( Pl Bg) - gB( Pl Bf))(x). 

On the other hand, let us consider a P-sequence with base point x. Then we can 
assume by a continuity argument that x t+ i x t-\ f° r ^ = ^i^ c -\- I n tfi a t case, 
since x — xq, we can rewrite the right hand side as follows 



RHS 



f(xo) r , s g(xi)-g{x ) , g(x )-g(x-i) 



Xl/2 - X-l/2 



Pi \ x l/2) Pl\F-\/2 



X\ — Xo Xq ~ X-i 



9{xo) r , v f(xi)-f{x ) , f(x )-f(x-i) 

Pl(Xl/2) pl(x_ 1/2 ) 



Xl/2 - X-1/2 



x\ — xq xq — X-i 



Since the terms in f(xa)g(xo) cancel, we easily obtain the desired relation. □ 
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Theorem 3.16. Let us consider a ^-sequence (xj) and a non-negative integer m. 
Then, if p is regular at xq,x\, . . . ,x m and ^1(^-1/2) = Pi{x m+ i) = 0, then the 
eigenfunctions fo, fx, ■ ■ ■ , f m satisfy the following orthogonality relations 



(fjjk) = fj( x i)fk(xi)p(xi)(x i+ i - Xl _i) = for j ^ k. 

i=0 

Proof. Lemma \'S. 151 gives 

tn 

(f,Lg)-(Lf,g) = £ D( Pl W(/, g))(x t )(x l+ 1 - x t _ 1 ) 



i=0 



i=0 



Pi(x)W(f,g){x) 



0. 



Taking f = fj, g — fk and recalling that /ij 7^ /i^ for j 7^ fc, the result follows. □ 

Remark 3.17. The inner product described in Theorem 13 . 1 61 can be complex. 

Remark 3.18. In Section 0] we shall study how free we are at the time of choosing 
the functions p and pj. Using this it can be checked that, when 

oui x _ 1/2 (pi) > ord xo Oi) = ord xl (pi) = ■ • • = ord Xm (pi) < ord Xm+1/2 (pi), 

we can modify the election of p and p\ so that hypotheses of Theorem 13 . 1 61 hold . 

4. Explicit formulas for the canonical forms 

In this section we apply our methods to study the hypergeometric polynomials 
that arise when P takes one of the following canonical forms: T, G, Q, A or C. 
Although the canonical forms O and E could also be treated with our methods, we 
shall not cover them here, as we explained at the beginning of Section [21 In each 
case we shall follow the following steps: 

1. Parameterization of a and r in terms of certain (Laurent) polynomial(s). 

2. Expressions for a.%, (3\, o~j, Tj and pu in terms of our parameterization. 

3. Taylor formulas for the eigenfunctions fo, fi, ■ ■ ■ fd- 

4. Rodrigues formulas for the eigenfunctions fo, fi, ■ ■ ■ fd- 

5. Classical families of hypergeometric polynomials which belong to that 
form. 

When dealing with Rodrigues formulas, we shall give in each case a method to 
construct the function p. The iterated functions pj can be easily obtained from 
p by applying Theorem 13.101 On the other hand, we shall assume again in what 
follows that the eigenvalues pk are pairwise distinct. 

4.1. Continuous case. V(x,y) = x 2 + y 2 — 2xy and 6t — xo- 

1. In this case we do not need any parameterization. 

2. It can be checked that we have 

Pk - m = (i - k)(a 2 (i + k) - a 2 + 0i) 
(Tj(x) = a(x), Tj(x) = t(x) + ja'(x). 
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3. Taking xq so that a(xo) = 0, we obtain from Theorem 13.41 

d 3 f k M=d kfk ( k )ii tm+^m 



j J Al a 2 (k + i) - a 2 + 0\ 
1 — j 



4. As we have pointed out in Remark |3.14l the Rodrigues formulas for this case 
are given by the solutions of the functional equations 

^ - (x) — p(x)r(x) and Pj(x) = p(x)a 3 (x). 

In general, we can not avoid to obtain ramified solutions. Quite surprisingly, this 
is the unique canonical form satisfying this property. Namely, as we shall see 
below, for any other canonical form we can always find non-ramified solutions of 
the corresponding functional equation for p. 

5. Taking <r(x) = x 2 — 1 (resp. <j{x) = x or u(x) — 1), we obtain the Jacobi (resp. 
Laguerre or Hermite) polynomials. Note that, if a(x) = 1 for all x £ C, there 
are no roots of a. In other words, the Hermite polynomials do not arise from the 
expressions given above and other methods are required to obtain them. However, 
the corresponding Rodrigues formula provided by Theorem 13 . 1 01 and Remark 13. 141 
also holds in this case. As we shall see in our analysis for the T canonical form, 
the g-Hermite polynomials behave in the same fashion. 

4.2. Arithmetic case. P(x, y) — (x — y + \){x — y — |) and 9 t = xq ± t. 

1. Let us consider two polynomials x ± ( x ) — 72^ 2 + l^x + with the same 
coefficient for x 2 . Clearly, there exists a bijective correspondence between these 
kind of pairs and P2 [x] x Pi [x] given by 

+ r x ( x , T ( x ) , x X + ( x ) + X~0) 
X^{x) = <j{x) + a{x) = 

X~{x)=a{x) — t(x) = x + (x) - X~{x)- 

2. Wehave/x fe -^ = (i - k)(j 2 (i + k) + 7^ ~7f - 72), a 2 =72 and ft =7f f -7f- 
Taking the P-function 9 t = t, we easily obtain the relations a^{t) = x ± {t)- Thus, 
it turns out by Lemma 13.21 that &f(t) — X ± (^ i §)■ Therefore, Oj and 73 are 
parameterized by 

X f{x)= X ± (x±l). 

In particular, we can write Tj(x) = x + ( x + §) — X~( x ~ §)" 

3. We have to consider two essentially different kind of P-functions: 

(a) If 9 t — X{) + t, we have (?~{t) — X~( x o + t). Thus, we take xq such 
that X~i x o) = to satisfy the hypothesis of Theorem 13.41 In that case, 
T i{6i/2) = X + ( x o + i) an d Theorem 13 . 41 gives 



jj fj-j 72 [i + k) + j{ - 7! - 72 
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(b) If 8t — xq — t, we have a (t) = x + i x o — t). Thus, if we choose the point 
xq so that X + ( x o) = 0, we obtain (6 , i / 2 ) = — X~{ x — i)< In particular, 



9jfk(0j/ 2 ) = d k fk 



k-l 

n 



-X Oo - i) 



72(i + fe) +7i -7i 



72 



Every case with pairwise distinct eigenvalues gives rise to an expression of this kind. 

4. Taking the P-function 6 t = t, the Rodrigues formulas given by Theorem 13.101 
arise from the solutions of 

p{t) X + {t) = p(f+l)x"(t + l) 

Non-vanishing entire solutions are available whenever are not identically zero. 
Moreover, those solutions can be obtained by multiplying the suitable solutions of 
the following root cases 



X+(x) = 1 
X~(x) = x-£~ 



x 0) 



x ■ 
1 



p(x) 



p(x) 



1 



T(l 



C" 



1 



x+(x) = c^o 

X~{x) = 1 
X +(x) = 1 

x -(x) = c^o 



p{x) = C 



p(x) = C x . 



r(i - x + £+) 

Other solutions are available multiplying p by a 1-periodic meromorphic function. 

5. Taking x^{ x ) — i x ~ ^t){ x ~ ^t)i we obtain the Hahn polynomials 

-k, £ - - x,k - 1 + f£ + £f - £+ - 

Co — Co ) Co — Cl 



-k, x £, 



o • 



1 



-.0 



(a) A(i)=A(^) 3 F 2 

(b) A(x) = / fc (#) 3 F 2 

4.3. Quadratic case. P(x, y) = x" + y _.. ., . h . 

1. Now we consider a polynomial x(i) = 74^ + 73^ + 72^ + 71 i + 70- Then, the 
parameterization is given by 

x(t) = a{t 2 )+tr{t 2 ) 



Co Co > Cl Co 
2 - 2a;y - \{x + y) + i and t = (i + i ) 2 - 



a(t 2 ) 



x(*) + x(-*) 



and 



xft) - x(-*) 

2f 



2. We have 012 = 74, A = 73 and 

Pk - Pi = (i - k) (74(1 + k) - 74 + 73) • 

Taking 8 t = t 2 and applying Lemma 13.21 we get 5^ (i) = x(±^ + §) f° r 3 > 0. 
Hence, <7j and Tj are parameterized by the polynomial 



Xj( x ) = xi x +h)- 



In particular, we have 



X (t+l)- x (-t+l) 



' 2t 
It suffices to consider the P-functions 9t - 



(t + to) 2 . In that case, we have 
&~{t) = x(—t " h) and so we impose x(~ *o) = 0. This condition provides the 
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factorization = (t + to)xo(t) for some xo G ^3[x]. Therefore, it turns out that 
Ti(0i/ 2 ) = Xo(ta + i). Theorem El gives 

\jj 74(« + k) - 74 + 73 
Every case with pairwise distinct eigenvalues gives rise to an expression of this kind. 
4. Taking the P-function 6t — t 2 and following Theorem 13. 101 we have to solve 
p(f) X (t) = p{(t + lf) X (-t-l) 

Non-vanishing entire solutions are always available. Namely, any solution arise as 
a product of functions like 

Pit 2 ) = ' 



r(i-£ + t)r(i-e-f)' 

which solves the functional equation for x(t) — t + £ and C. Moreover, other 
solutions appear when we multiply p by any even 1-periodic meromorphic function. 

3 

5. Taking x( x ) = ]^[0 T + £u) an d ^o = Coj we obtain the Wilson polynomials 



f k (t 2 ) = h(e ) 4 f 3 



Co + t, £o - t, *-l + fo + fi + 6 + Cs 

& +6, Co + 6 



4.4. Some remarks on a functional equation. Finally, it remains to study the 
geometric and trigonometric canonical forms. In order to give explicit Rodrigues 
formulas for these particular cases, we shall need to solve functional equations of 
the following type 

(9) F(x) = R(x)F(qx), 

where R stands for a non-zero rational function. An open subset £1 of the complex 
plane will be called g-invariant when qfl = O. We are interested in non-zero 
meromorphic solutions defined in a g-invariant connected open set CI. We shall 
need to apply the following basic results. 

Lemma 4.1. Let Fj and F 2 be non-zero meromorphic solutions of the functional 
equation Q in a q-invariant domain CI with respect to the rational functions Ri 
and R 2 respectively. Then F]T 2 and F1/F2 are non-zero meromorphic solutions of 
@ in CI with respect to the rational functions R1R2 and R1/R2 respectively. 

Lemma 4.2. A function F solves ij^jl in a q-invariant domain CI for some q G C 
if and only if F solves the functional equation 

F(x) = S(x)F(q- 1 x) in Cl with S(x) = -. — - . 

R{q~ L x) 

Lemma 4.3. Let us suppose that a given q-invariant connected open set Cl contains 
the point 0. Then, if there exist a non-zero meromorphic solution of @ in Cl, it is 
unique up to a constant factor. 

Proof. The quotient F of any two meromorphic solutions of the functional equation 
(PI satisfies F(x) — F(qx) in Cl. Therefore, it suffices to write this relation in terms 
of the Laurent series of F in a neighborhood of 0. This completes the proof. □ 
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By Lemma [4.21 the case \q\ > 1 can be reduced to the case \q\ < 1. Now we 
provide some particular solutions for \q\ < 1 which will be useful in the general 
case. Namely, the g-shifted factorials 

g^(x) — (£;x; q)oo with H(x) = 1 — £x and CI = C 
h^(x) = (q/£x; q)^ with ~R(x) = — and fi = C \ {0}. 

Given a non-zero rational function R, we can always write R in the form 

with r € Z, ( a non-zero complex number and £k,Vk G ^ for 1 < fc < n . Then, 
using the particular solutions and h^, Lemma f4. II assures that the function 

g-yh-, 



gihi 



(gihiY n^/fe with 7 = (-i) r C 
/c=i 

solves the corresponding functional equation. F is a meromorphic function in C\{0}. 
In fact, if £ = 1 and r — 0, F is meromorphic in C and we know by Lemma 14 . 31 that 
F is the unique meromorphic solution in C up to a constant factor. 

Remark 4.4. Let G be meromorphic in C and biperiodic with periods 2m and 
logg. The function F defined by F(expy) = G(y) is meromorphic in C \ {0} 
and solves the functional equation F(a;) = F(qx). The theory of elliptic functions 
provides infinitely many functions of this kind. In particular, this shows there is no 
uniqueness of meromorphic solutions in C \ {0} of the functional equation (jSJ). 

n 

Remark 4.5. By the previous discussion, when q n = J^[(6c/?7fc) the function 

fe=i 

is elliptic with periods 2wi and logq. Moreover, E has a zero at — log^ and has a 
pole at — log77fc for 1 < k < n. The rest of zeros and poles of E can be obtained by 
periodicity. Therefore we have obtained, in our terminology, a proof of a classical 
existence result. This explicit expression allows to explore the consequences derived 
from the non-uniqueness of solutions p and pj of the Rodrigues formulas. We shall 
give more details on this topic in the following Paragraphs 4.5 and 4.6. 

Remark 4.6. The functional equation Q has non- vanishing rational solutions for 
q = e 2 W™ whenever the following extra hypothesis hold 

fc=0 

These solutions could be used to obtain Rodrigues formulas for such values of q. 
However, we have decided not to include them in this paper for lack of space. 

Remark 4.7. Putting f(t) = F(q t ), © takes the form f(t) = R(<?*)/(t + 1). This 
kind of equations have been solved in for \q\ < 1 in terms of the r g -function. 
In general, these solutions can not be regarded as non-ramified solutions in q l . 
This alternative functional equation is also considered in where some solutions 
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f(t) = F(q t ) are given. However, the authors seem not to be specially interested on 
solutions which are functions of the variable q l and they do not observe the role of 
elliptic functions here. For more on the connection between g-calculus and elliptic 
function theory, we refer the reader to the references £Q and [Hj. 

4.5. Geometric case. P(a;, y) — x 2 + y 2 — 2axy and 9 t = q^XQ. 

1. Let us consider two polynomials x ± ( x ) = l2 x2 + 7i x + 7o with the same value 
at x = 0. Then, the parameterization is given by 

X {%) = o-(x) + xt(x) a(x) = 

— / \ f s A -A" 1 x + (x) - X~(x) 

X {x)=a{x) xt(x) t(x)= ( X _ x -ij x ■ 

2. We have a 2 = ( 7 + + 7 j)/2, ft = ( 7 + - 7 2 ~)/(A - A" 1 ) and 

\i— k \k~i 

^ - f M = {x _ x r } 2 (7 2 + A t+fc - 1 - 7 2 -A 1 - l - fc ). 

Taking 9 t =q t , it turns out that cF ± (t) = x ± (l t ) an d of{t) = X ± (^ ±i Q t )- In other 
words, the polynomials Oj and Tj are parameterized by 

X f(x)= X ± (^x). 

In particular, we obtain 

X +(A^)-x'(A-^) 
T ^ ) = (A-A^ ■ 

3. We have to consider two essentially different kind of P-functions: 

(a) If 9 t = q l XQ, we have &~{t) = x~(q tx o)- Thus, we require x~( x o) = 
to satisfy the hypothesis of Theorem 13.41 Under this assumption, we have 
Ti(0i/2) = X + ('? I a ; o)/(A — X~ 1 )X l Xo and Theorem 13. 41 gives 



d j f k {e j/2 ) = d k f k [ k .] J -^)-^l[ 



fc-i 



x + (q tx o)/ x o 
i " 72" -i2 ( i 1 ~ i ~ k 



(b) If 9 t = q t xo, we have a (t) = x + (q t %o)- Thus, we take xo such that 
X + ( x o) = 0. This gives Ti(6>i/ 2 ) = -A 1 x _ (g _i a:o)/(A - A -1 ):^- In particu- 
lar, 



d j fk(6 j/ 2)=d k f k \ k ] g-[(5)-a)] J] 



k-1 



J 



-x (q tx o)/ x o 

1 X=j 72~ ~ l2l 1 ^^ k 



Every case with pairwise distinct eigenvalues gives rise to an expression of this kind. 
When xo = 0, the terms x ± (q ±Zx o)/ x o must be replaced by the obvious limits. 

4. Taking the P-function 9t = q* and following Theorem I3.1UI we have to obtain 
non-zero meromorphic solutions of 

p(x)x + (x) = p(qx)x~(qx) 

p J+1 (x) = p J (\^x) X ± (\ ±U - 1) x). 
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in a P-invariant domain f2. Non-vanishing meromorphic solutions in C \ {0} are 
available whenever \q\ < 1 and are n °t identically zero. It suffices to multiply 
the solutions of the following root cases 



where m — m_ — m + and 7 = ( — l) m (l~ q m ~ / C, + . Besides, Remark 14.41 provides 
infinitely many meromorphic solutions in C \ {0}. For particular cases, choosing 
the simplest solution might not be completely trivial. 

5. Taking X ± (x) = (1 - $x)(l - tfx) with £± ^ and ^ 0, we obtain the 
q-Hahn polynomials 



(a) fk(x) = f k (l/Za) 302 ' 9 W/$>6 



(b) /*(*) = /*(!/#) 302 ( ^ ^ «* 



4.6. Trigonometric case. P(x,y) = x 2 +y 2 —2axy+a 2 — l and 6 t = ^ U ° ^ — ?fo__ 

1. In this case, the Laurent polynomial x(u) = 7-2"~ 2 + 7-iu~ 1 + 70 + 71 u + 72^ 2 
provides the following parameterization 

, . fu + u~ 1 \ (A — A _1 )(u — u~ v ) [u + u~ l \ 
= °{—2-) + " 4 2~ ) 

, u + u -^ = X (u)+x(u-i) ^ r /" + «-^ = rt")-x(Q 
V 2 / 2 V 2 / (A - \- r )(u - u- 1 ) 

2. We have a 2 = 2(-f 2 + 7-2), Pi = 4(72 - 7-2)/(A - A -1 ) and 

Mfc - Mi - 4 A ( ' A ^ A 1 " )2 t ( 72 A''+ fc - 1 - 7 - 2 A 1 - l - fe ). 

Taking 9 t = {q* + g~*)/2, we obtain 5^(4) = x{q ±l ) and 3^(i) = x(A J g ±4 ). That 
is, the polynomials Oj and Tj are parameterized by 

Xj(u) = x(A 3 m). 

In particular, we have 

(A — A _1 )(u — u^ 1 ) ' 

3. It suffices to consider the P-functions 8 t = (g'uo + q~ tu o 1 )/'2- In that case, we 
have &~(t) = x(? - * u o )• Thus, we take Uq so that xi^ 1 ) = 0- This condition 
provides the factorization x(u) — (u—Uq 1 )u~ 2 xo(u) forsomexo € Ps[a;]. Therefore, 
it turns out that 



t. 



2X0(9^0) 



X 3l u 2 (X- A- 1 )' 
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Theorem 13.41 gives 

djfk(0j/ 2 ) = d k fk 



9 11 2(72 - 7-29 1 - 1 " 

? — j 



Let us notice that, for x( u ) = 7±2M ±2 , our methods provide a Rodrigues formula. 
However there is not any uq such that x( u o ) = 0- Hence we can not apply Theorem 
13.41 These cases correspond to the g-Hermite polynomials. 

4. Now we apply the ideas of Paragraph 14.41 to discuss some Rodrigues formulas 
for the trigonometric case when |q| < 1. We begin by recalling that the relations 

n=|— : ueCl'j and Cl' = |tt : — ^ 6fi] 

provide a bijective correspondence between the class of P-invariant open sets Cl of 
the complex plane and the class of A-invariant open subsets Cl' of C \ {0} satisfying 
that u £ Cl' if and only if u^ 1 6 Cl'. Moreover, given Ci and CI' as above, the relation 

establishes a bijection between meromorphic functions G in Ct and meromorphic 
functions F in CI' satisfying F(it) = F(u _1 ). Consequently, we shall work with 
functions pj of the form 



i ( U+ 2 ) = M M ) 



Pj I 

where hj is meromorphic in CI' and hj(u) = hj{u^ r ) for all u G fi'. Taking the 
P-function 9 t = (q* + q~ t )/2 with u = q l , we have to solve 

h(u)x(u) = M^M? -1 " -1 ) 

We can always write the Laurent polynomial \ in the form 

r— 1 

x(«) = C«-]I( 1 -^ u ) 

t/=0 

where r — s = max{fc : 7^ ^ 0} and £ = 7fc such that A;o = min{fc : 7^ 7^ 0}. 
Using the arguments of Paragraph 14.41 and assuming \q\ < 1, we observe that the 
functions 

h i {u) = W^qU(i^qU 
are meromorphic solutions in C \ {0} for x( u ) = 1 — £u an< i invariant under the 
change u 1— > On the other hand, when x( u ) = the system above reduces 
to 

h(u) = qu 2 h(qu) 
h j+ i(u) = A 1 ^M ±1 /i J (A ±1 w). 

Let h ± (u) = (tA 1 / 2 ^; A) 00 (tA 1/2 /w;A) 00 . We have h ± (u) = ±\ 1 ' 2 uh ± {\u). In 
particular, the functions h + and h~ solve the functional equation h(u) = qu 2 h(qu). 
Moreover, if we define recursively the constants k^ by 

kf +l = ±X^ j kf with = 1, 
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it can be easily checked that the functions = k^h^ are meromorphic solutions 
in C \ {0} for \( u ) — an d invariant under the change u i— > u . In summary, 
by Lemma |4.1I we have meromorphic solutions in C \ {0} of the system above for 
any Laurent polynomial \ 

r — 1 r — 1 

hj(u) = (kfh±y(u) n hf{u) = {±iy s \^-^)\h ± y(u) J] hf(u). 

i/=0 u=0 

Therefore, the functions pj obtained by this process are meromorphic in the complex 
plane. For instance, if we have s = 2 and r = 4, we can choose the solution h^h~ 
for x( u ) — u an d so we get 

(io) Pj ( u + U ~ 1 ) = (-A)yG) . ( Xu2 '^ x / u2 ^U . 

As in Remark 14 .41 we can find other meromorphic solutions in C \ {0} of the given 
functional equation. Namely, it suffices to multiply ft by a function / satisfying 
f(expx) = g{x) with g even, 27ri-periodic, log g-periodic and meromorphic in C. 

Remark 4.8. The solutions we use for the functional equation h(u) — qu 2 h(qu) do 
not follow the general solution provided in Paragraph l4.4l The main motivation for 
our election has been the simplicity of the resulting functions pk ■ 

5. If 7_2 and 72 are non-zero, we obtain the Askey- Wilson polynomials. Namely, 
normalizing if necessary, we can write 

3 

X (u) =u- 2 l[(l-^u) 
in such a way that 7-2 = 1 and 72 = £o£i£2£;3- Therefore we get 



= fk 



Co + Co 1 A ^ f q k , Cow, Zo/u, Co6669 fc 1 
V C0C1, C0C2, C0C3 



5. Concluding remarks 



q, q 



1. In the canonical forms analyzed in Section the eigenfunctions are families 
of polynomials depending of five parameters which correspond to the coefficients 
of a and r. The number of free parameters can be reduced by considering two 
proportional hypergeometric operators to be equivalent. Moreover, recalling how L 
is transformed under the action of Aff(C) 

g LV 1 = C^ 9 (D sP ) 2 + Cr 9 S 9P Df p , 

we can use the isotropy subgroups of the „4i-orbits provided in Table I to identify 
different hypergeometric operators by suitable affine transformations. In particular, 
we can reduce the number of essential parameters in each canonical form. Namely, 
we have two free parameters in the continuous canonical form, three parameters 
in the arithmetic and geometric cases and four parameters in the quadratic and 
trigonometric canonical forms. 

2. The Rodrigues formula given in [2] for the Askey- Wilson polynomials arise from 
(|TUj) by taking 



pf w (u) = k(u)p 3 (- 
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Here, the function k is defined as follows 

w , = (qu 2 ;q) 00 (q/u 2 ; q)oo(u - v," 1 ) 
(Xu 2 ;q) 00 (X/u 2 ;q) 

DC 

The function k is meromorphic in C \ {0} and satisfies the functional equation 
k{qu) — k{u). However, we have k{u) = —k(u~ l ). Therefore, the functions used in 
[2] are ramified in the variable 

x=-(u + u ). 

The function k can be easily regarded in terms of Jacobi's sn elliptic function. 

3. Let us consider a function e meromorphic in C \ {0} and such that 

o The function e satisfies e(u) — e(qu) — —e(l/u). 

o The function e{\u)p\{\{u + it -1 )) is regular for X < \u\ < A -1 . 
Then it can be proved that, assuming that g, S (0,1) and adapting the arguments 
of |141 3.10.5.2], the hypergeometric operator L is symmetric with respect to the 
product 

(f,g) = J f(x)g(x)p{x)e(x)dx 

where e(cos(s)) = e(e ts ) for < s < ir. This gives the orthogonality relations for 
the Askey- Wilson polynomials fk- The function k defined on the previous point 
satisfies the required conditions on e. In fact, by the Ricmann-Roch theorem, it can 
be checked the it is unique up to a constant factor. Hence, this allows to recover 
the well-known orthogonality relations given in |2] . 
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